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ABSTRACT 

We infer the normahzation and the radial and angular distributions of the number density of satel- 
lites of massive galaxies {logiQ[AI^/MQi] > 10.5) between redshifts 0.1 and 0.8 as a function of host 
stellar mass, redshift, morphology and satellite luminosity. Exploiting the depth and resolution of the 
COSMOS HST images, we detect satellites up to eight magnitudes fainter than the host galaxies and 
as close as 0.3 (1.4) arcseconds (kpc). Describing the number density profile of satellite galaxies to be 
a projected power law such that P{R) oc K^^, wc find 7p = —1.1 ± 0.3. Wc find no dependency of 7p 
on host stellar mass, redshift, morphology or satellite luminosity. Satellites of early-type hosts have 
angular distributions that are more fiattened than the host light profile and are aligned with its major 
axis. No significant average alignment is detected for satellites of late-type hosts. The number of 
satellites within a fixed magnitude contrast from a host galaxy is dependent on its stellar mass, with 
more massive galaxies hosting significantly more satellites. Furthermore, high-mass late-type hosts 
have significantly fewer satellites than early-type galaxies of the same stellar mass, possibly indicating 
that they reside in more massive halos. No significant evolution in the number of satellites per host 
is detected. The cumulative luminosity function of satellites is qualitatively in good agreement with 
that predicted using subhalo abundance matching techniques. However, there are significant residual 
discrepancies in the absolute normalization, suggesting that properties other than the host galaxy 
luminosity or stellar mass determine the number of satellites. 



Subject headings: galaxies: dwarf - galaxies: evolution - galaxies: formation - 
- gravitational lensing 
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1. INTRODUCTION 

Cold Dark Matter (CDM) simulations of structure for- 
mation successfully match observations of both the spa- 
tial and mass distribution of super-galaxy scale struc- 
ture. However, they are less successful at matching 
observations of structure at smaller scales. For exam- 
ple, these simulations over-predict the number of low 
mass sate llite com panions of Milky Way mass halos 
dKauffma nn et al.l[! 993: Klv pin et aLlll99a iMoore et ahl 
[19991: iStrigari et al.ll2007) . This discrepancy may demon- 
strate a fundamental breakdown of the ACDM paradigm 
at small scales, possibly indicating that assumptions 
about the nature of the CDM particle are incorrect. For 
instance, DM particl es may have high e r kinetic energy 
than expected (e.g. iCoh'n et all l2000t iSchneider et al.l 
120111 : iMenci et al.H2012[ ) leading to suppressed fluctua- 
tions in the primordial power spectrum at small scales. 
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there by dec reasing the numbe r of present day satellites 
(e.g. iKamionkowski fc Liddld [2000; Zentner fc Bullock 
120031) . Alternatively, the missing satellite problem may 
be due to the difficulty of observing very faint satellites 
due to low star formation efficiency in low mass halos. 
A variety of baryonic processes can produce suppressed 
star formation in lo w mass halos (e.g iTliovi l fc Wcinbcr 
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1996; Gncdin 2000 
2010; .SpringeL201G[ ). However, it is not easy to obser- 
vationally confirm which processes have the strongest ef- 
fects, in part due to the small number of faint satellite 
galaxies studied outside of the Local Group. 

In the Local Group, the luminosity function of satellite 
galaxies has been stu died to a faintness of My = —2 (see 
iKoposov et al.l 120081 and references therein) . Outside 
of the Local Group, measurements have been restricted 
to brighter satellites. Multiple studies using data from 
SDSS have measured the lu minosity f unction of satel- 
lites at low redshift (iGuo et al. 201 lal iLiu et al.l 120111 : 
iLares et al.ll20TlL iStrigarTfc WechsleJl2011| ). The deep- 
est study by IGuo et al.l (|2011al) . measured the luminos- 
ity function reliably for sate llites up to 7.5 m agnitudes 
fainter than the hosts, while IStrigari fc Wechsl cr (201l]) 
placed upper limits on satellite numbers for satellites up 
to 10 magnitudes fainter. All of these studies agreed on 
two main conclusions. The first was that the luminosity 
function of satellites is dependent on the host luminos- 
ity. This is due to the fact that for massive host galax- 
ies, the host stellar mass is believed to be a non-linear 
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function of the host virial mass Ce.g.lBehroozi et al.|[20lol: 
iWake et al.ll2011l:lLeauthaud et al.i201lD . Thus^ile the 
number of satelhtes at a fixed virial m ass ratio with th e 
host is expected to be scale invariant (|Kravtsovll201Cl[ ). 
the luminosity ratios may not be. The second conclusion 
from these studies was that while the number of faint 
satellites appears to be consistent with that in the Lo- 
cal Group down to an absolute V band magnitude of -14 
(IGuo et al.ll2011aHLares et al.ll2011 l :lStrigari fc Wechsleii 
1201 on the bright end only 10-20 percent of hosts with 
Milky Way and Andromeda luminosities host satellites 
as luminous as the Large and Small Magellanic Clouds 
(L/SMCs). The rarity of high mass companions is consis- 
tent with predictions from ACDM simulations, indicating 
that the Milky Way luminosity function is an outlier on 
the bright end rather than that there is an issue with the 
theory. However, a significant discrepancy between the- 
oretical predictions and observation of the satellite lumi- 
nosity function only begins to appear for satellites with 
absolute magnitudes fainter than My > —14. Thus, in 
order to test whether the missing satellite problem is due 
to the Milky Way having an anomalous lack of satellites 
on the faint end, it is necessary to study even fainter 
satellites. 

In addition to providing a key to understanding star 
formation in low mass halos, and potentially allowing 
us to constrain the mass of the DM particle, satel- 
lite galaxies are believed to have played a significant 
role in the evolution of the size of the most massive 
galaxies (M* > 1O^°M0) (Bovlan-Kolchin ct al. 2006; 
IBovlan-Kolchin fc Mai [20071 : [Hopkins et al.ll20TfP . The 
role of minor mergers in the evolution in size and 
mass has been observed in numerous observational 
studies which used pair counting and assumptions about 
merger time-scales to estimate the m inor merger rate 
(iLe Fevre et al.|[200C 
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20081: [Bundy et al.l 
Robaina et al.[ [201C 
201 It ). Furthermore, measurements of disturbed mor- 
phology and color gradients in an elliptical hosts 
as evidence for re cent mergers ([Kavirai et al.l [2009[ : 
[Kavirai et al.[[20ll[) show that minor merging is a key 
contributor to low level star formation seen in the 
outskirts of early-type galaxies. These studies have 
all been limited to the study of the most massive 
companions with stellar masses at least ten percent of 
their host mass, which is ten times higher than the mass 
ratio between the Milky Way and the LMC. In order 
to better constrain the effect of minor mergers on the 
evolution of massive galaxies, it is necessary to push the 
study of companions at cosmological distances to lower 
stellar masses Csee also [Bundv et all [2007[ : [Naab et all 
[200l iFakhouri et al.[[2010ll 

The spatial distribution of satellites about their host 
galaxies is intimately tied to the luminosity function 
of those satellite galaxies. Two competing effects con- 
tribute to the link between satellite luminosity function 
and spatial distribution. On the one hand, baryons are 
thought to play an important role in the preservation 
of subhalos in the inner regions of the host dark matter 
halo by steepening the total mass profiles and making 

* But see also lNipoti et al.l ||20T1 ) 



them less susceptible to tidal stripping. Numerous hy- 
drodynamical and semi-analytic simulations have found 
that subhaloes with baryons are more centrally located 
in the host dark matt er halo than pure dark matter 
subhalos (e.g. Blumcnthal et al.l 119 861: iLibeskind et al 
20101 Romano-Diaz ct al. 2009; We inberg et al.l [200 ' 
Maccio et al...2006: .Gao et al.,.2004) . 



At the same time, baryons will also steepen the mass 
profile of the host, increasing the strength of tidal 
shock heating of the gas and stars in satellites as they 
pass near the center of the halo, with the amount 
of heating and mass l oss strongly dependent on the 
orbits of the satellites (iGnedin et a,l.l [19991: [Choi et all 



200l [Bodan-Kolchin fc Mal l2007l: IDolae et al.l 1200911 
D'Onghia et all (120101 ) find that the effects of these in- 
teractions in suppressing star formation is stronger for 
less massive subhalos. For these reasons, the slope of 
the radial profile of satellites is of key importance to any 
analysis attempting to recover the physics that governs 
the interactions between satellites and their host halo 
and central galaxy. 

The radial profile of the number density of satellites 
has been measured with somewhat contradictory results. 
Assuming a single power-law model for the projected 
number density of satellites as a function ofdis- 
tance from the host, where N{R) cx RJ, iChenI (|2008[ ) 
measured 7^ — —1.5 ± 0.07 for satellites brighter than 

Mr 17.5 in SDSS. In a study of SDSS satellites 

brighter than Mg < -2 1.2, W atson ct al. (2010) mea- 
sured 7^ = -1.2 ± 0.1. In [Nierenberg et all (|201lD (here- 
after Nil), we measured the slope of higher redshift satel- 
lites up to 5.5 magnitudes fainter than their hosts (on 
average about My < —16) to be 7^ = — l.O^oi- 

There are several likely explanations for the discrepan- 
cies in these measurements. One possibility is differ- 
ences in the technique used to separate satellites from 
background/foreground interlopers. Alternatively the 
difference in inferred slopes may be due to differences 
in host masses studied, or differe nces in the luminos - 
ity of the satellites. For example, I Watson et al.l ([20111 ) 
measured the radial profile of satellite galaxies brighter 
than Mr < — 18 between 7 and 280 kpc as a function 
of satellite luminosity, and found that the faintest satel- 
lites had significantly shallo wer radial profiles than their 
more luminous counterparts. iTal et aL[(|2012[ ) also found 
fewer faint satellites in groups within 25 kpc for satellites 
brighter than Mg < —22.8. However outside of this re- 
gion, they found that all satellites followed the same ra- 
dial profile which was well described by a combined Ser- 
sic-|-NFW pr ofile, mimicking the total mass pro file which 
follows (jGavazzi et al.l 120071 : [Auger et al llMol) for 
massive galaxies. In contrast to these studies. iGuo et all 
(2012) find that for satellites brighter than rur < 22 in 
SDSS, fainter satellit es are more centrally co ncentrated 
than bright satellites. iBudzvnski et al.l (|2012|) also found 
an excess of fainter satellites in the innermost regions for 
satellites in groups and clusters. 

To make progress on understanding the satellite radial 
profile, it is necessary to make measurements that take 
into account systematics such as the host mass and satel- 
lite luminosity. In addition, in order to learn about 
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trends in merging rates i t is interesting to study evolu- 
tion in the radial profile. iBudzvnski et al.l ()2012D found 
no significant evolution in the radial profile of satellites 
in groups and clusters between redshifts 0.15 and 0.4. 
In this work we extend the measurement of the satel- 
lite radial profile to fainter satellites and to higher red- 
shifts. 

The angular distribution of satellites also contains infor- 
mation about the effects of anisotropic, filamentary ac- 
cretion in ACDM and the shape of the host dark matter 
halo. Pure dark matter simulations predict that satel- 
lites should be found in an anisotropic distribution which 
is aligne d with the major axis of the host dark matter 
halo (e.g. iZentner et al. 2005; Zcntncr 2006; Kncbc et al. 
[200llAubert et al.ll2004: ,Fahenbacher et al...2007. ,2008. ) . 
Thus by studying the orientation of satellite galaxies 
with respect to the major axis of the host light profile 
it may be possible to learn about the relative orienta- 
tion between the stars in the host galaxy and the dark 
matter halo they reside in. In the Milky Way and An- 
dromeda, satellite galaxies appear in an anisotropic dis- 
tribution, prefe rentially aligned with the minor axis of 
their hosts (e.g. lMetz et al.ll2009[ ). When averaged, how- 
ever, studies of satellites in SDSS find that satellit es are 
randomly located around la te-type galaxies (Baili n et al.l 
I2008t lYegorova et~alll2011[ ) and aligned with the major 
axis o f the host light profile for early -type hosts (jBrainerdl 
120051 lAgustsson fc BrainerdI [20101 ). The studies of the 
satellite angular distribution outside of the Local Group 
have all been limited to relatively bright satellites (within 
1-2 magnitudes of the host magnitude), and they have 
not attempted to measure the relative ellipticity of the 
satellite angular distribution. As stated above, measur- 
ing the angular distribution of fainter satellites is of key 
interest because several simulations predict that satel- 
lites that have very elongated orbits will lose the most 
gas as they pass through the host galaxy. 

Finally, the measurement of the luminosity function and 
radial profile of satellites at intermediate redshifts can 
be combined with gravitational lensing stu dies to con- 
strain star form ation efficiency in subhalos (jTreul [2010l : 
iKravtsovl I2010L and references therein). Recently low 
mass subhalos have been det ected in mass reconstruc- 
tions of gravit ational lenses (Vegctti et aL| 120121 120101 : 
iMcKean et"all 2005; MacLeod ct al,, 200^ By analyz- 
ing the probability of d etecting satellites of these masses 
near the lensed images, iVegetti et al.l (|2012| ) constrained 
the low mass end of the mass function of subhalos to have 
a slope of a = l.l^g ^. At present, the combination of 
gravitational lensing analyses with luminosity function 
measurements is the only method to directly measure 
stellar to virial mass ratios for satellites outside of the 
Local Group. 

In this paper we address three questions: 1) What is 
the spatial distribution and cumulative luminosity func- 
tion of the satellites of massive galaxies? 2) How does 
the spatial distribution vary with satellite luminosity? 
3) How do these properties vary with host morphology, 
stellar mass and redshift? We conclude by discussing 
how our results can be interpreted in the context of 
ACDM. 

In Nil, we developed a method of modeling and sub- 



tracting the host light profile to allow us to detect 
faint {rrisat — 'nihoat > 5-5) satellites at intermediate 
redsh ifts near early-ty pe host galaxies in the GOODS 
(Giavalisco et al."2004') fields. We also introduced a sta- 
tistical model which we used to simultaneously infer the 
spatial and angular distribution of satellites along with 
their numbers. In this work, we expand our analysis to 
include the satellites of host ga laxies selected from the 
COSMOS ijScoviUe et al.l |2007[) field which is approxi- 
mately 20 times larger than the GOODS field although 
somewhat shallower. This dramatic increase in our sam- 
ple size allows us to analyze the properties of the radial 
and angular distribution and the number of satellites, for 
satellites more than a thousand times fainter than their 
host galaxies. Furthermore, the greatly increased num- 
ber of host-satellite systems allows us to analyze these 
properties in bins of redshift, host morphology, satellite 
luminosity and host mass. This is of key importance be- 
cause, as has been discussed, the number and luminosity 
of satellites are believed to depend strongly on the prop- 
erties of host galaxies. In order to make our results more 
easily relatable to theoretical predictions, we select our 
host samples by stellar mass and luminosity and analyze 
the number density of the satellite radial profile in units 
of i?2oo in addition to units of the half- width at half-max 
of the host light profile. 

This paper is organized as follows: In fj2l we describe 
the images and catalogs used in this analysis. In ij3]we 
discuss the selection and properties of the host galaxy 
sample. In 2] we briefly review the host modeling and 
subtraction technique developed in Nil. In iJ5]we present 
the binned radial and angular profiles of satellites to pro- 
vide a qualitative sense of the signal. In [J6] we review 
the methodology for a full Bayesian statistical analysis 
of the data developed in Nil which allows us to use a 
catalog of object positions and magnitudes to infer the 
presence of satellites against the homogeneous signal of 
background/foreground objects. In fj7|we discuss the the- 
oretical model of the cumulative number of satellites per 
host a s a function of host luminosity from iBusha et al.l 
(poll which we compare with our measurements. In 
Sj8]we explain how we tested our statistical analysis and 
its limitations. In §[5] we present the results. In [TUl we 
compare our results with similar studies performed at low 
redshift. In § [Tl] we discuss the broader implications of 
our results. In §[12] we provide a concise summary. The 
Appendix contains more detailed explanations of many 
of the methods used in this paper. 

Throughout this paper, we assume a flat ACDM cos- 
mology with h = 0.7 and Q.-m = 0.3. All magnitudes 
are given in the AB system (jOkdllQ?^ unless otherwise 
stated. 

2. IMAGING AND CATALOGS 

To achieve the depth and area required to study satellite 
properties and their evolution, we use imaging andpho- 
tometric catalogs taken from the COSMOS survey [j. In 
Nil we used the deeper GOODS survey to develop and 
test our method. However, the GOODS survey is not 

^ COSMOS photometric catalogs and surveys are available at 
http: / /irsa.ipac.caltech.cdu/data/COSMOS/datasets.html 
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sufficiently wide to allow us to perform a binned statis- 
tical analysis of the satellite population as a function of 
host properties as we do in this paper. 

The COSMOS field has extensive spectroscopic and pho- 
tometric ground-based coverage which we use to cat- 
alog the redshift, stellar mass and morphology of the 
host galaxi e s. Spe ctroscopic redshift measurements from 
iLillv et all ()2007f ) are available for roughly 1/7 of the 
host galaxies. When spectroscopic measurements are not 
available we use ground-based photometric redsh ift mea- 
surements from the catalog bv lllbert et al.l (|2009( ). which 
is the same catalog we use for all stellar mass measure- 
me nts. Finally for morphological categorization, we use 
the iCassata et al. (2007) morphological catalog. The 
iCassata et al.l (|2007) catalog is created by a computer 
algorithm which is known to fail in certain cases. For 
instance, spiral galaxies with large bulges can be mis- 
classified as ellipticals. Thus we visually confirmed all 
morphological classifications using COSMOS ACS imag- 
ing. This resulted in the re-categorization of ~ 10% of 
the morphologies of host galaxies. 

3. HOST GALAXY SELECTION 

In selecting our host galaxy sample, we had two main 
goals. The first was to ensure that we selected hosts that 
maximized our ability to detect neighboring satellites, 
and the second was to select a host sample that could be 
easily related to predictions from simulations. In Nil, 
we achieved the first goal by selecting massive early-type 
host galaxies which had relatively smooth light profiles 
that made it easy to find nearby companions. However, 
galaxy morphology is n ot well reproduced by current sim- 
ulations (e.g. Hoyle et al. 2011), making it difficult to 
select matching host galaxies from simulations. In con- 
trast, simulations have had success matching the stellar 
mass function o f massive galaxies, to dark matter halos 
in simulations (iBerlind fc Weinberg l2002t iZheng et al.l 
[20071 IConrov fc Wechsledl2009n . making a mass-selected 
host sample more appealing. 

To compromise between the two goals, we selected rela- 
tively high stellar mass \og[M^/MQ\ > 10.5 host galax- 
ies, as these have well studied stellar masses, redshifts 
and luminosities, and a greater number of sat ellites vis- 
ible down to a fixe d app arent magnitude (e.g. iKravtsovl 
I2010t iBusha et all 120 111 ). We study hosts in a redshift 
range between 0.1 < z < 0.8. This range allows us 
to study evolutionary properties while still guarantee- 
ing that we can detect satellites with luminosity contrast 
from the hosts equivalent to that of the SMC relative to 
the Milky Way at all redshifts. By doing this, we expect 
to observe approximately one satellite per host. We ex- 
clude z < 0.1 galaxies, which are few and too extended 
in angular size to analyze in the same way as the more 
distant sample. 

We also required host galaxies to be relatively isolated 
in order to ensure that they themselves are not satel- 
lites of larger central galaxies. This is important for our 
analysis because we do not want to count objects associ- 
ated with the larger host as satellites of a smaller galaxy. 
Using the stellar mass and redshift catalogs, we include 
only host galaxies which are not within the i?200 of a 
neighbor that has more than its stellar mass and is at 



the same redshift within measurement uncertainties. We 
calculate R20 Q given stellar mass using Equation 3 from 
iDutton et al.l ()2010[ ). which provides a by eye fit to the 
observed and inferred relationships between stellar and 
halo mass of galaxies as a function of morphology. For 
early-type galaxies the best fit function is: 
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And for late-types: 

» = io"(to^)"'[5 + 5(to^^)J 

Where y = (M200) /M* and x = M*. 

Although there is uncertainty in the stellar mass esti- 
mate and additionally in R200, we consider these values 
to be know with absolute precision for the purposes of 
our analysis. We tested the effect of uncertain stellar 
mass on our analysis by doubling masses before calculat- 
ing i?2oo and found no significant impact on our inference 
result. The final sample has 1901 early-type and 1524 
late-type galaxies. The distribution of stellar masses, 
redshift and absolute r-band magnitudes are shown in 
Figure [TJ 

As discussed in the Introduction, there is theoretical and 
observational evidence that satellite properties may de- 
pend strongly on the properties of host galaxies. In order 
to study these trends, we divide host galaxies into bins 
of redshift, morphology and stellar mass when perform- 
ing our analysis of the satellite population. We choose 
two bins of stellar mass with 10.5 < logj^Q < 11.0 and 
11.0 < logigM^ < 11.5, and two bins in redshift with 
0.1 <z< 0.4 and 0.4 <z< 0.8. 



When we analyze the satellite population, we study bins 
of cumulative magnitude contrast from the host galaxy, 
such that all satellites have magnitudes brighter than 
'm-sa.t — "T-host < Am. In each bin of Am we only study 
host-satellite systems such that the host is at least Am 
magnitudes brighter than the limiting survey magnitude. 
This means a different subset of host galaxies is used to 
study Am = 8.0 satellites than Am = 2.0 satellites, lead- 
ing to some variation in host properties within a stellar 
mass, redshift and morphology bin. In Tables [2] and |3l 
we summarize the average stellar mass, redshift and lu- 
minosity of the host galaxies used in each Am analysis, 
split into bins of redshift, stellar mass and host morphol- 
ogy. Host galaxies used to study the faintest satellites 
tend to be at lower redshift than the average in a fixed 
redshift bin. 

4. DETECTION AND PHOTOMETRY OF CLOSE 
NEIGHBORS 

In order to study the satellite population, we require an 
accurate catalog of object positions and magnitudes as 
close as possible to the host galaxies. Companions of 
bright galaxies are difficult to study because they are in- 
trinsically faint and often obscured by the host galaxy 
light. This is a serious issue when attempting to mea- 
sure the slope of the power law-radial profile of the satel- 
lite spatial distribution, as the innermost regions of the 
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Fig. 1. — The distribution of stellar mass, redshift and absolute r band magnitude for upper: early and lower: late-type hosts. 
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system provide the best constraint on the slope. Fur- 
thermore, close to the host, the ratio of the number of 
satellites to background/foreground galaxies is the most 
favorable. In Nil, we developed a method of removing 
the smooth component of the host galaxy light profile to 
overcome some of these challenges. This process makes 
automated object detection much more accurate and re- 
liable near the host galaxy. We use the results to update 
the COSMOS photometric catalog with newly detected 
objects, and to replace the photometry for objects that 
had already been detected near the hosts. 

In this section we briefly review the method of host light 
modelling and subtraction developed in Nil and we dis- 
cuss the new objects this method allows us to identify. 
All SExTRACTOR parameters used are listed in Table S] 
in the Appendix. 

4.1. Companion masking and host light subtraction 

Host subtraction allows us to accurately identify faint 
objec ts near our ho st galaxies using similar SExtrac- 
TOR ([Bertin &: Arnouts ,1996 ) parameters to those that 
were to make the COSMOS photometric catalogs. These 
parameters were conservative as they were intended to 
study objects across a large range of angular sizes and 
morphologies. To ensure that we are able to identify faint 
objects that the conservative parameters might have oth- 
erwise missed, we remove the light from the images in a 
three step process. 

In the first step, we use SExtractor 
(|Bertin fc Arnouti 119961 ) to identify and then mask 
objects in a small cutout region of 10 i?h near the host 
galaxy, where i?h is the second moment of the host 
galaxy light profile along its major axis (AJMAGE). 
In this step, we select more 'liberal' SExtractor 
parameters than are used in the final round of object 
detection, erring on the side of masking noise, to ensure 
faint satellites are masked. Two separate rounds of 
object identification are combined to make the final 
mask, one which is tuned to detect diffuse objects farther 
from the central galaxy and the second which is tuned 
to detect faint compact objects closer to the central 
galaxy. Finally, we use a two dimensional, elliptical 
b-spline model to model the masked image and subtract 
the host light pr ofile. The model we use is similar to 
that described bv lBolton et al.l (|2005l[2006l) . 

4.2. Object detection and photometry 

We detect objects in the host-subtracted images using 
SExtractor parameters tuned to match the object de- 
tection and photometry that was used to make the COS- 
MOS catalogs. Appendix [B] contains a comparison of our 
photometry and object detection to that of the COSMOS 
catalogs in a large non host-subtracted field. 

In order to test our sensitivity to low surface brightness 
objects, we simulated faint sources near our hosts with 
Sersic indices and effective radii and tested our recov- 
ery rate and photometric accuracy for these objects af- 
ter host subtraction. Results of these simulations can 
be found in Figure [TT] in the Appendix. We use the re- 
sults from these simulations to identify a minimum radius 
at which we can detect at least 90 percent of simulated 



satellites with accurate photometry, and define this as 
the minimum radius at which we study the properties of 
the satellite population. We find that we can accurately 
recover satellites with MAG_AUTO Isu < 25.0 as close 
as 2.5 i?h in COSMOS for the majority of hosts with 
early-type light profiles. This corresponds to a mean dis- 
tance of 1'.'2 (7 kpc) with a standard deviation of 0'.'7 (3 
kpc). The inner detection boundary is slightly higher for 
host galaxies with > 11.0 and z < 0.4 as these tend 
to have light profiles that extend above the background 
further from the host centers. Late- type galaxies have 
more extended light distributions in addition to spiral 
arms which are difficult to distinguish from neighbor- 
ing galaxies, thus we choose a more conservative inner 
boundary of 4 _Rh(3±l arcseconds, 17± 6 kpc) for these 
hosts. On average, these minimum radii correspond to 
0.02 and 0.07 i?2oo for early and late-type galaxies re- 
spectively. 

4.3. Properties of Objects Detected in Cutout Regions 

In this section we compare the properties of newly de- 
tected objects after host light subtraction to the proper- 
ties of objects already in the COSMOS photometric cat- 
alog near the host galaxies. The upper panel of Figure 
[2] shows the distribution in the contrast in MAG_AUTO 
measurement {Sm= m — rrihH between hosts and ob- 
jects detected within 2'.'5 of the host galaxies. The dis- 
tributions are compared for objects already in the COS- 
MOS catalog and newly detected objects. For both types 
of objects, photometry is performed after host light re- 
moval. Newly detected objects are about a magnitude 
fainter than objects that were already in the COSMOS 
photometric catalog, with average i5m values of 3.22 com- 
pared to 4.44 for previously detected objects within the 
same region, with typical measurement uncertainty of 
0.05 mag. 

The lower panel of Figure [2] shows the number density 
of objects as a function of distance from the hosts. The 
number density of objects in the COSMOS photomet- 
ric catalogs drops within the inner 1'.'5, while the num- 
ber density of newly detected objects rises, more than 
doubling the number of COSMOS detected objects in 
this region. The sum of the two number density sig- 
nals increases steadily with decreasing distance from the 
host galaxy. Thus host subtraction and rigorous mea- 
surements of completeness in the innermost regions are 
necessary for an accurate measurement of the radial pro- 
file and number of objects close in projection to the host 
galaxies. 

5. FIRST LOOK 

Before describing our model for the radial and angular 
profiles of objects near the host galaxies, it is instruc- 
tive to show these distributions in spatial bins in order 
to provide a visual representation of the data. How- 
ever, binning is inherently limited because it requires the 
averaging of data, thereby losing information. Further- 
more it is not conducive to accurate subtraction of fore- 

Note the use of lower-case <5m, which denotes a specific con- 
trast from the host and is different from Am which describes the 
allowed maximum contrast between host and neighboring objects 
for a particular data set. 
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distance scaling; the number density of sources increases 
as a power-law near the hosts. At large radii, the number 
density becomes dominated by the isotropic and homo- 
geneous distribution of objects not associated with the 
hosts, represented by the gray dashed lines. 

In Section [6] we describe how we analyze the number 
density signal by inferring the combined properties of 
the sat ellite and background/foreground populations. In 
Section there is a comparison of the results using the 
two distance scalings. 
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Fig. 2. — Comparison of the properties of objects detected in the 
COSMOS catalogs to those of newly detected objects after host 
light subtraction. Upper: The distribution of magnitude differ- 
ences from hosts {5m= m — m^) within 2'.'5. Lower: The number 
density of objects as a function of distance from hosts. Newly 
detected objects are closer to the hosts than those in the COS- 
MOS photometric catalog and make a significant contribution to 
the measurement of object number density within 2 arcseconds 
from the center of the host. 

ground/background galaxies over a range of redshifts. 
Thus we do not perform our analysis on the spatially 
binned data, but instead use this section to justify our 
model choices in Section [01 

5.1. Distance Scaling and Radial Distribution 

We scale measured object distances to account for the 
range of redshift and host mass scales in our sample. 
A scale relating to the host light profile is the natural 
choice for the observer as this will vary with host redshift 
as w ell as host ma. ss according t o the siz e-mass relation 
(e.g. iTruiillo et a ll [2006 : Williams et al.| [2010). For this 
distance scale, we use i?h which is AWIN_IMAGE from 
SExTRACTOR . We also perform a parallel study with all 
distances scaled by i?200 of the host galaxies. Unlike i?h, 
i?200 can be calculated in dark matter only simulations 
and is thus a better choice when attempting to compare 
results with simulations. However, estimating i?200 re- 
quires multiband photometry and stellar mass modeling 
which is not always possible, so it is useful to perform this 
analysis using both distance scalings to see if one choice 
or the other leads to systematic differences. 

Figure[3]shows the average number density of objects as a 
function of distance from the hosts, with distances scaled 
by i?h in the upper panel and i?200 in the lower panel. 
The behavior is qualitatively similar for both choices of 




100 




R/R,„ 



Fig. 3. — The average number density of objects near hosts as 
a function of radial distance. Upper: In units of the second order 
moment of the intensity of the host-light profile along its major 
axis (i?h): and lower: In units of R200 estimated from the host 
stellar mass. 



5.2. Angular Distribution 

In Figure 2] we show the angular distribution of objects 
within 10 i?h, plotted for all hosts, early-type hosts and 
late-type hosts, where = is aligned with to the major 
axes of the host light profiles. This figure only includes 
host galaxies with axis ratio b/a<0.6, to ensure that the 
direction of the host major axis is clearly measurable. 
As background/foreground objects are expected to be 
distributed isotropically relative to the host galaxy, any 
anisotropy we observe is caused by correlated structure 
presumably, in the form of satellites. In this region we 
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expect a significant contribution to the number density to 
come from satellites, as evidenced by the strong satellite 
signal within this region in the upper panel of Figure 

13 

Near early- type galaxies, the angular distribution shows 
a dominant component aligned with 9 = 0. A 
Komogorov-Smirnoff KS test gives a probability of 
10^* that the objects near early-type galaxies have a uni- 
form angular distribution. In contrast, the objects near 
late-type galaxies appear more isotropically with a KS 
probability of being uniform of a few percent. This sim- 
ple examination has been done without any effort to re- 
move background/foreground contamination. However, 
it indicates that the satellites of early and late-type hosts 
may have different angular distributions relative to their 
host light profiles. To test this further, we separate satel- 
lite populations based on host morphology in our infer- 
ence of the parameters of the satellite spatial distribution 
in Section [6lTl 



6. 



JOINT MODELING OF SATELLITE AND BACKGROUND 
GALAXY POPULATIONS 



In this section we review how we infer the number of 
satellites and their radial and angular distributions from 
a catalog of object positions and magnitudes. The ba- 
sic procedure is to model the observed positions and 
number density of objects as a combined signal from 
background/foreground objects and satellite galaxies. 
These two populations have significantly different spa- 
tial distributions; satellite galaxies increase in number 
near the host galaxies and for early-type galaxies ap- 
pear preferentially aligned with the major axis, whereas 
background/foreground objects have a homogeneous and 
isotropic number density signal. By inferring the proper- 
ties of the combined signal, and using prior information 
about the background/foreground objects, we isolate the 
satellite number density signal. 



In Subsection 16 . 1 1 we discuss the details of this model for 
the satellite spatial distribution and in Subsection 16.21 
for the background/foreground objects. Theoretical and 
observational justifications for the model choices, along 
with a schematic illustrating a possible realization of our 
model can be found in Nil. 

6.1. Satellites 

6.1.1. Satellite Spatial Dtstnhution 

We construct our spatial distribution model to investi- 
gate three main components. The first is the slope of the 
radial profile of the satellite number density. We model 
the projected satellite radial distribution as a power-law 
P{R) oc . The second two parameters describe the 
ellipticity and orientation of the angular distribution of 
satellites relative to the host light profile. We model the 
satellite angular distribution as elliptical with axis ratio 
Qs between its major and minor axes and with an orien- 
tation, \(j)\ of its major axis relative to the major axis of 
the host light profile |(/)| is allowed to vary between and 
90 degrees, where indicates alignment with the major 
axis of the host light profile. 

Given these parameters, the probability of finding a satel- 
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Fig. 4. — The average number of objects at angle \9\ from the 
major axis of the host Hght profile, within 10 -Rh, for hosts with 
elongation (b/a) less than 0.6. The distribution of objects near 
early-type hosts is more aligned with the host light profile than 
the distribution of objects near late-type hosts. 

lite at a position (i?, 0) is given by: 



K<^[cos^ {e-(j)) + l/ql sin^ (61 - (j))X'^'^RdRd6 



(3) 



The normalization of the angular part of this distribution 
is given by a generalized elliptical integral. 

For some of the host galaxies such as face-on spirals and 
very symmetric early-types, the orientation of the host 
major axis is not well defined. To account for this in 
our inference, we infer the relative fiattening between 
satellites and host galaxies. 

Defining the ellipticity of the satellite population to 
be: 



Cosmic evolution of substructure 



9 



We relate the ellipticity of the host light profile to the 
ellipticity of the satellite angular distribution via the pa- 
rameter A: 

= 1 + 

The parameter A can be understood as follows: when 
A is zero, satellites are always distributed isotropically 
around host galaxies, regardless of how flattened the host 
light profiles are. When A is one, the satellite population 
follows an ellipse with the same elongation as the central 
light profiles, as A goes to infinity, the distribution ap- 
proaches a straight line. In this work, we only examine 
whether the satellite angular distribution is rounder or 
flatter than the host light profile, so we restrict our in- 
ference to values of A between zero and two. 

6.1.2. Number of Satellites per Host 

To ensure completeness, we study the satellite population 
brighter than a fixed magnitude contrast from the host 
galaxies mobj — < Am, where the host galaxies in a 
Am bin are all at least Am brighter than the background. 
We model the number of satellites Ng in a given Am bin, 
between iJmin < -Rsat < ^^max SiS constaut for all hosts in 
that bin. 

Both i?inin and i?max are chosen based on observational 
considerations. i?min is the limit where accurate satel- 
lite detection is possible given host light subtraction, as 
described in Section S) i?max is the limit where a satel- 
lite signal is apparent above the background in Figure [3l 
which is 45 Rh for early and late- type galaxies. These 
outer limits correspond on average to 0.5 and 1.0 i?200 
for early and late- type galaxies. We use these average 
values for i?max when performing the inference in units 
of -R200 ■ 

6.2. Background/ Foreground Objects 

We model the background/foreground number density 
signal around each host as a homogeneous, isotropic 
signal, with mean surface density Eh with magnitudes 
brighter than the field magnitude limits (/814 < 25.0). 
We inform our inference of E;, by measuring the local 
background around each of our hosts and using the mean 
and standard deviation of this distribution to create a 
prior. 



We follow the method recommended by 'Che n et al.l 
(|2b06. ) of estimating the local background around our 
host galaxies (rather than taking the field average) in 
order to accurately measure E^. This is important for 
removing correlated line of sight structure which is not 
within the virial radii of the host galaxies. In Nil, we 
measured the background between 45 and 70 R^, where 
no significant satellite signal was apparent. For consis- 
tency, in this analysis, we estimate the background in 
the same region when scaling distances by i?h- When 
performing the analysis in units of i?200j we calculate 
the background betwee n 1.0 and 1.5 ji!2 nn, following the 
simulation results from lLiu et al.l ()201l( l. 



mean and standard deviations separately for each data 
set that we run the inference on and use this local back- 
ground information as a prior when inferring the prop- 
erties of the satellite distribution. 

For a given value of Am and host galaxy magni- 
tude, mh, only a fraction of the total number of back- 
ground/foreground objects will appear in our analy- 
sis. To calculate this fraction, we model the cumu- 
lative distribution function (CDF) of the backgroun d 
number counts by a power- law fe.g. 'Beni tez et al.ll2004f) . 
Defining Eb,o as the background density measured for 
all objects with magnitudes brighter than the limiting 
survey magnitude mii,„, the number density of back- 
ground/foreground objects around the j'*'' host is given 
by: 



(6) 



We measure the background slope to be 0.305±0.005 near 
low mass early and late-type galaxies and 0.300 ± 0.005 
near high mass host galaxies. 

7. THEORETICAL PREDICTION FOR THE NUMBER OF 
SATELLITES PER HOST 

Our choice to model the number of satellites as being 
constant for all hosts within a fixed magnitude contrast 
from the host magnitude is a simplification of a more 
complicated picture. Dark matter only simulations pre- 
dict that dark matter halos should host an approximately 
constant number of dark matter subhalos with a given 
dark matter mass fracti on of the host halo, regardless 
of the host halo mass (jKravtsovl 120101 and references 
therein). This scale- invariance relates in a non-trivial 
way to our chosen observable which is the number of 
satellites we expect to see within a fixed magnitude con- 
trast from from host halo. This is due to the non-linearity 
of the stellar-mass to halo-mass relationship for hosts 
with M* > 1O^°M0 as can be seen, for example, in 
iBehroozi et all (|2010l ). 



Various studies have shown that the local projected num- 
ber density of objects is correlated with host stellar mass 
and morphology. To account for this, we measure the 



In recent work, iBusha et al.l ()2011[ ) (hereafter Bll) 
modeled the number of satellites per Am m bin as 
a function of host luminosity using SubHalo Abun- 
danc e Matching (SHAM) techniques (jKravtsov et al.l 
120041 : IBehroozi et al.l 120101) . to c onnect dark-matter 
halos in the Bolshoi simu lations, (jKlvpin et al.l 120111 : 
iTruiillo-Gomez et al.l 1201 ID to the r-band luminosity 
function of galaxies. In order to properly mimic the ob- 
servational selection function in this work, we have re- 
produced the measurements of Bll using the luminosity 
function from the AGES simulations (Cool et al 2012) 
as applied to the Bolshoi halo catalog at the appropriate 
redshift. 

To test how the inference method responded to the non- 
Poissonian distribution of satellites predicted by theory 
we ran simulations using the distribution of host lumi- 
nosities in our sample to generate the number of satellites 
per host based on Equation 8 from Bll. Spatial po- 
sitions for satellites and background/foreground objects 
were drawn stochastically from the model described in 
Section ini As desired, our inference accurately returned 
the mean number of satellites per host along with the 
other input parameters in the simulations. 
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TABLE 1 
Summary of model parameters 



Parameter 


Description 


Prior 


Satellite Model 






Ns 


Number of satellites per host 


U(0,20) 


Tp 


Logarithmic slope of the satellite radial distribution 


U(-10,0) 


A 


Relative flattening between satellite angular dist. and host light 


U(0,2) 


l</'l 


Angle between the major axis of satellite angular dist. and that of the host light. 


U(0,7r/2) 


Background Model 


Number density of all background objects with 7gi4 < 25 
Logarithmic slope of the background number counts 


varies with host mass and Am 
N(0.300,0.005)/N(0. 305,0. 005) 



b 



U(a,b) denotes a uniform distribution between a and b. 

N(/i,cr) denotes a normal distribution with mean /i and standard deviation a. 



8. ANALYSIS 

In Nil, we provide details for the construction of the 
posterior probabihty distribution function (PDF) which 
ahows us to infer values for the parameters of our model 
given the data and prior knowledge of the background 
listed in Table [TJ For each parameter in each subset of 
the data listed in Tables [5] and [31 we compute the poste- 
rior PDF using a Markov Chain Monte Carlo (MCMC) 
method. At least 10^ iterations per chain are performed 
in order to ensure convergence. 

To study variation with host properties, satellites are 
analyzed in bins of 'high' and 'low' host stellar mass 
corresponding to 10.5 < \og^Q[M^/MQ\ < 11.0, 11.0 < 
log;^o [^^h /^^0] < 11-5) as weU as 'low' and 'high' host 
redshift corresponding to 0.1 < z < 0.4 and 0.4 < z < 0.8 
and early and late- type host galaxies. Every bin in Am 
is analyzed separately for each bin in host morphology, 
redshift and stellar mass, using an appropriate prior on 
Sb,o estimated using the local background for that subset 
of host galaxies. 

In order to combine results from different data sets, we 
bin the posterior PDFs for the parameters of interest and 
multiply them together. When data sets were analyzed 
in different regions around the host, we first re-normalize 
the inferred satellite numbers using the posterior median 
values of 7p and Ng to account for the differences in ex- 
amined areas, before combining the posteriors. 

We tested our inference by running simulations with 
varying background and satellite properties. We find 
that in order to accurately infer the radial profile 7p we 
need a minimum of 50 host galaxies and at least 20 satel- 
lites. The first requirement guarantees that the inference 
has an adequate estimate of the background density, and 
the second that there is a sufficient satellite signal. Infer- 
ring the angular distribution of satellites is more difficult 
as this requires two dimensions of information. We find 
that at least 50 satellites are necessary to recover A and 
Furthermore, the inference on \(j}\ becomes inaccu- 
rate for values of A less than 0.5, thus we do not report 
confidence intervals of |0| where the posterior favors val- 
ues of A less than 0.5. 

9. RESULTS 

We divide this section into the three characteristics of the 
satellite population that our inference constrains, namely 
the radial and angular distribution of satellites as well as 



the cumulative luminosity function. For each character- 
istic we discuss variation with host morphology redshift 
and stellar mass as well as with satellite luminosity by 
comparing results in high and low redshift and stellar 
mass bins as well as in bins of host morphology, as dis- 
cussed in Section [8l Summaries of the inference results 
for each subsample Tables [2] and |3] in the Appendix. 
To show typical results. Appendix |D] contains the full 
bivariate posterior PDFs for all model parameters for 
high mass, early and late-type hosts for the Am — 4.0 
bin. 

9.1. Radial Distribution 

In Figure [5] we show the inferred value of the projected 
slope of the radial profile of the satellite number density 
as a function of Am, divided by host morphology, stel- 
lar mass and redshift. We further compare these results 
when distances are scaled by i?h and by i?2oo- 

This figure contains several important results. The first 
is that there is a significant detection of a population of 
objects with a power-law radial distribution up to 6.5 
magnitudes fainter than low redshift hosts, and 5.5 mag- 
nitudes fainter than high redshift hosts. Using the aver- 
age luminosity of the host galaxies in each of those data 
sets, this corresponds to approximate absolute r band 
magnitudes of -16.1 and -17.4, respectively. The former 
is about one magnitude brighter than Sagittarius, the 
latter is similar to the absolute magnitude of the present 
day SMC. 

Second, there is no significant difference in the inference 
of 7p when distances are scaled by i?h or by R2oo- This is 
remarkable given the number of assumptions that go into 
estimating the virial mass given the stellar mass, and is 
a very useful result for comparisons between theoretical 
and observational work for future satellite surveys where 
accurate host stellar masses are not available. For sim- 
plicity, for the remainder of the paper we discuss results 
only for distances scaled by i?2oo- Results for both dis- 
tance scalings are listed in Tables are listed in Tables [5] 
and [3] in the Appendix. 

Third, there is no significant variation in 7p as Am varies 
with fixed host properties over ranges as large as 5.5 mag- 
nitudes. This is consistent with the physical processes 
governing the satellite luminosity function being not spa- 
tially related to the host galaxy over this fairly large lu- 
minosity range within the distance scales we study, given 
our measurement uncertainties. 



Cosmic evolution of substructure 



11 



The average values of 7p for host galaxies of all morpholo- 
gies are: -l.liO.l, -1.3±0.4, -1.2±0.1 and -l.OiO.l 
for satellites of low mass-low redshift hosts, high mass- 
low redshift hosts, low mass-high redshift hosts and high 
mass-high redshift hosts respectively. 

Individual subsets of host morphology do not show any 
significant deviation from these valuesEj As there are no 
significant differences in the results for these samples, we 
argue that the satellite population can be well described 
by a power law with 7p — —1.1 ± 0.3 within the level of 
precision afforded by our data. 

9.2. Angular Distribution 

The satellites of early and late-type galaxies display 
markedly different angular distributions where the sam- 
ple of satellites and host galaxies is large enough to allow 
an inference on the parameters on the angular distribu- 
tion. As an example, in Figure [5] we show the two dimen- 
sional bivariate posterior distributions for the parameters 
A and |(/)| for satellites with Am < 4 for early and late- 
type hosts with 0.4 < z < 0.8 and 11.0 <Log[M*]< 11.5. 
The results for this subsample which are representative 
of all non-uniform posterior PDFs. One sided confidence 
intervals for A and \4>\ for other non-uniform subsamples 
are listed in Tables[2]and[3] Satellites of early-type galax- 
ies show strong anisotropy, with the most likely value of 
A at 2 and A > 0.9 at a 68% confidence level. This in- 
dicates that the angular distribution of satellites is more 
flattened than the light profiles of the early-type galax- 
ies. Furthermore, the satellites of early-type galaxies, are 
found preferentially along the major axis of the host light 
profiles. The offset, \4>\ has a most likely value of and 
is less than 20 degrees at a 68% confidence level. 

In contrast, the average angular distribution of the 
satellites of late-type galaxies is consistent with being 
isotropic relative to the host disks, with A having a most 
likely value of and A < 1.1 at a 68% confidence level. 
As discussed in Section [HI the inference on \(f>\ is unreli- 
able for small values of A as |(/)| has decreasing influence 
on the likelihood function for small values of A. 

There is a concern with regards to describing the elliptic- 
ity of the angular distribution of satellites as a fraction, 
A, of the host ellipticity (see Equation [5]). Namely, late- 
type galaxies are much more elongated than early-type 
galaxies but there is no reason to expect the satellites 
to reflect this. It may be that the satellites of late- 
type galaxies have anisotropic distributions on average 
but that the flattening of their distribution is less than 
the extreme flattening of the disks, making them look 
isotropic in comparison. To test whether the inferred 
isotropy is due to the extreme flattening of late-type 
hosts, we re-ran the inference on late-type hosts, this 
time artiflcially changing any late- type axis ratio b/a to 
0.6 that was previously flatter than 0.6. We found that 
the posterior PDF still favored A = 0. However, the in- 
ference on the angular distribution for late-type hosts is 
difficult due to the large region obscured by spiral arms. 

Although the high mass low redshift sample of late- type galax- 
ies appears to have a steeper satellite number density radial profile, 
this sample is relatively small (only 54 hosts) and does not shovf a 
statistically significant deviation. 



To try to maximize the signal to noise, we performed a 
combined inference on late- type hosts of all redshifts and 
stellar masses. In this case, there was evidence for slight 
anisotropy, with satellites aligned with the galactic disk. 
Thus, an analysis of the satellites of late- type hosts with 
a large data set is warranted before any strong conclu- 
sions are reached. 

9.3. Cumulative Luminosity Function 

We detect a significant population of satellites as faint 
as 6.5 magnitudes fainter than their host galaxies. Re- 
sults for the inferred number of satellites per host as a 
function of maximum contrast between host and satel- 
lite magnitude (mhost — Wsat < A^) can be found in 
Tables [2] and [3] A dominant uncertainty in this analysis 
is caused by the covariance between 7p and Ns as can be 
seen in Figure fT2l the same total number of objects can 
be achieved if there are more satellites following a shal- 
lower radial profile or if there are fewer satellites with a 
very steep radial profile, and more background objects. 
We alleviate this degeneracy by performing the inference 
a second time, using a Gaussian prior on 7p from the re- 
sults in Section 19.11 Given that 7p showed no significant 
variation given measurement uncertainties as a function 
of satellite magnitude, host redshift or host morphology 
(see Section [^?T|) we apply the same Gaussian prior on 7p 
to all data sets, with mean —1.1 and standard deviation 
of 0.3. 

Results for the inferred number of satellites using the 
prior are consistent with results without the prior and 
are listed in Tables [2] and [3] The prior on jp allows for 
a significant detection of satellites up to 8 magnitudes, 
or more than a thousand times fainter than their host 
galaxies for the low reshift host sample. Using the aver- 
age absolute r band magnitude of hosts, this corresponds 
to satellites with Mr ~ —13.5/ — 14.7 for low and high 
mass host samples respectively, and is similar to the ab- 
solute magnitude of Fornax. The range of Am does not 
change for the higher redshift sample but the average 
measurement uncertainty is decreased. 

In Figure [7] we plot the number of satellites per host, Ns, 
in increasing bins of Am, for varying host morphologies, 
stellar masses and redshifts between 0.07 and 1.0 i?200 
(^ 17 — 200 kpc). The points for 'AH' morphology hosts 
come from binning and multipl ying the posterior PDFs 
of the early and late-type hosts E3- We always combine 
the posteriors even when the inference on the late-type 
host galaxies is not sufficiently determined to plot (i.e. 
for Am > 4 for more massive, high redshift late- type 
hosts). The satellite CLF is a fairly constant power-law 
with slope approximately N{L) oc 7j-0-5±o.i^ There is 
a slight upward shift in the normalization for satellites 
with Am > 6.5, which is attributable in part to the slight 
increase in average host stellar masses for these bins (see 
Tables [2] and [3]). We leave a more detailed unbinned 
analysis of the satellite luminosity function to a future 
paper. 

There is a strong dependence between early-type host 

This is valid as these data sets are independent. The product 
of the two posterior PDFs can be viewed as a measurement of the 
number of satellites per host for the entire sample 
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Fig. 5. — Projected radial profiles of the satellite number density distribution with distances scaled by Rjj and -R200 for satellites divided 
by host morphology, redshift and stellar mass. Horizontal lines indicate the average median value of 7p. Large markers indicate results for 
satellites of higher stellar mass hosts and small markers for the satellites of lower stellar mass hosts. 
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Fig. 6. — Bivariatc posterior distributions of |<^| and A for the 
satellites of Upper: early and Lower: late-type host galaxies, for 
hosts with 0.4 < 2 < 0.8 and 11.0 < log[A/*^^j/Af0] < 11.5 

stellar mass and number of satellites within a fixed value 
of Am. The dependence on stellar mass is less apparent 
for late-type host galaxies, with the more massive hosts 
having barely more satellites on average than less massive 
host galaxies despite the fact that the low and high mass 
late-type hosts have the same average stellar masses as 
the low and high mass early-type hosts. While there is 
a large difference between the number of satellites of the 
more massive early and late-type hosts, the numbers are 
the same for less massive hosts. 

Figure [7] also shows the Bll theoretical prediction for 
the mean number of satellites given the distribution of 
host luminosities in each sub-sample. The dotted line 
represents an extrapolation of the theoretical model, to 
luminosities where the simulation used begins to suffer 
from incompleteness effects. The theoretical model is 
in excellent agreement with observation for objects with 
0.1 < z < 0.4, while it tends to underpredict the abun- 
dance of satellites of higher redshift objects (although 
the slopes remains in good agreement). It should be 
noted, however, that the luminosity function used in the 
model in the range 0.4 < z < 0.8 suffers from significant 
observational uncertainties. Such uncertainties directly 
impact halo occupation in non-trivial ways, since param- 
eters such as and have a significant impart on the 
mass-luminosity relation for massive objects. Addition- 
ally, we assumed a fixed scatter in the mass-luminosity 
relation of 0.16 dex, which has been shown to be in good 
agreement with z = 0.1 galaxies (Behroozi et al 2010), 
but has yet to be explored at higher redshift. A full ex- 
ploration of the theoretical uncertainties relating to this 



prediction is beyond the scope of the current paper and 
is left for future studies. 

As expected, subhalo abundance matching predictions 
which are constructed without taking host morphology 
into account cannot capture the morphological depen- 
dencies. A possible explanation for th e variations is 
given by the morphology-density rela tion (jDressleiillQSOl : 
iPostman fc Gelleii llQSl iTreu et al.]|2003[ ). in the sense 
that at fixed luminosity, early-types tend to reside in 
denser environments and have more massive dark mat- 
ter halos than late- types. This in turn might affect the 
properties of luminous satellites and points to limitations 
of the present SHAM models. 

10. COMPARISON WITH PREVIOUS WORK 

We can make several comparisons between our work and 
low-redshift studies of the satellite population. To start, 
our inference on 7p — — 1.1 ± 0.3 is consistent with 
results from the low redshift SDSS study of LRGs by 
IWatson et al.l (|2010D and of satellit es with lumino sities 
similar to that of the LMC (ToUerud et al.ll2011[ ). It 
is slightly shallower, even th ough m arginally consistent, 
than what was measured bv iChenI ((200S) . Given uncer- 
tainties in the scale radii of host galaxies, we are unable 
to determine wheth er our resu l t is co nsistent with NFW 
as was observed bv lGuo et al.l ()2012D . Our result is con- 
sistent with our previous work (Nil), in which we stud- 
ied the satellites of early-type galaxies from the GOODS 
fields. 



In contrast to IWatson et al.l (|201lD . we do not find a sig- 
nificant trend between Am an d 7r< over the 5 magn itudes 
in Am that we studied. While IWatson et all (|2011l ) mea- 
sured the satellite radial profile changing from ^ —2 
for satellites with Mr = —21 to 7p ^ —1 for satellites 
with Mr = —18.5. The closest matching data-set in this 
study is the low-redshift, high-mass sample, for which 
we infer 7p ~ ^O.T^Q g for Am = 1.5 satellites (which is 
approximately Mr ^ —21) and — l.T^Qg for Am = 4.5 
satellites (approximately Mr ^ —18.5). The difference 
may be due to the fact that our enhanced detection tech- 
nique allows us to detect a significant new population of 
faint satellites which would have otherwise remained ob- 
scured, causing the radial profile to appear to be fiatter. 
Given our measurement uncertainties we cannot rule out 
a small increase in the concentration o f fainter sa t ellites 
towards the center as was observed by iTal et al.l (poTTI ) 
and lGuo et all (|20T1 . 

The dependence of the satellite angular distribution on 
host morphology is consistent with numerous previous 
low-redshift studies of bright satellites ( Bailin et aLll2008l : 
IBrainerdI 120051 : lAgustsson fc Brainerdll2010t) . which also 
found that satellites of early-type hosts tend to be found 
along the major axis of the host light profiles. Further- 
more, these works found no significant anisotropy in the 
angular distribution for the s atellites of late- type hosts 
(see also lYegorova et aLll2011|) . 

The satellite numbers we measure for lower redshift 
hosts are consistent with the results from low red- 
shift SDSS studies of hosts with similar masses. In 
Figure |51 we compare the number of satellites of low 
mass, low redshift host galaxies to the number in the 
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Fig. 7. — The cumulative number of satellites per host between 0.07 and 1 -R200 ^ function of the magnitude contrast between host and 
satellite galaxies, plotted for different samples of host redshift, stellar mass and morphology. Purple squares, red circles and blue diamonds 
represent all, early and late- type hosts respectively. The gray solid and dashed lines arc theoretical predictions for the satellites of high and 
low mass host galaxy samples respectively. Thin dotted lines indicate an extrapolation of the theoretical prediction which was made for 
satellites brighter than Mr < — 17. Note that the mean host stellar mass and reclshift within each bin shifts slightly towards higher masses 
and lower redshifts starting at Am = 6 (see Tables [2] and UJl . . 
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Milky Way and near o t her low mass, low r edshif t host 
galaxies by iGuo et alj (1201 lal) . iLiu et all (|2011[ ). and 
iStrigari fc Wechsleii (|2011[ ). The Milky Way has an ab- 
solute magnitude of Mr ~ —21.2, which is slightly lower 
than the typical host luminosity in the low stellar mass, 
low redshift data set in this sample {Mr ~ —21.6). How- 
ever we expect some passive evolution in the stellar mass 
to li ght ratio of order one between redshifts 0.8 and 0.1 
fe.g. ' Treu et al.|[2005[) . which would make the average 
luminosity approximately -21.3 at present day. Further- 
more, the typical stellar mass in the low mass, low red- 
shift subset of galaxies is log^g [M^/M©] ~ 10.7 which is 
equivalent to the stellar ma ss of the Milky Way within 
measurement uncertainties ()McMillan|[20ll . Thus the 
hosts are approximate Milky Way analogs at redshift 
0.3. 

In Figure [H we also plot two theoretical predictions. 
The Bll subhalo abundance matching is identical to that 
plotted in the corresponding panel in Figure [71 The sec- 
ond shows the 10-90 p ercent tails of the p redicted satel- 
lite distribution from iGuo et all ()2011bl) semi-analytic 
models of galaxy formation appli ed to the Millennium 
and Millennium II simul ations (jSpringel et al.l 120051 : 
iBovlan-Kolchin et al.ll2009h . 

At all comparable magnitudes our inferred number of 
satellit es i s consisten t with that me asured bvlGuo et al. 
2011ar ). ILiu et al.l (|20Tll) and IStrigari fc Wechsler 
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At the bright end, there are on average fewer 
satellites per host than what is observed in the Milky 
Way. The infrequent presence of Magellanic Cloud- 
equivalent satellites (Am = 2.0) of Milky Way mass 
hosts has been noted in numerous surveys before this 
(|Guo e t al. 2011a; Lares et al. 2011; Liu et al.ll20Tll ). We 
are able to measure the satellite number down to 8 mag- 
nitudes fainter than the host galaxies and corresponds to 
an average r band luminosity of about -13.3 which is sim- 
ilar to the Fornax satellite. At the faint magnitudes the 
satellite numbers fa l l just under the upper limits from 
IStrigari fc W echsler (2011) and the slope of the luminos- 
ity function is consistent with that at brighter satellite 
magnitudes. 



At higher redshifts, iNewman et al.l (|2011[ ) measured the 
pair fraction of galaxies Am < 2.5 between redshifts 0.4 
and 2. Scaled to the same region as our study, they found 
iVg = 0.48 ± 0.09. This number is somewhat higher than 
our measurement for low-mass host galaxies which had 
0.14 ±0.1 in the same region. This discrepancy is likely 
due to the fact that the authors used global rather than 
local background subtraction, meaning they estimated 
the average background across the entire field. This leads 
to a lower estimate of the number density of background 
objects than local background estimation as it does not 
take into account the filamentary clustering of galaxies. 
The authors estimate that using local rather than global 
background estimation increases their pair fraction by 
a factor of two relative to local background techniques, 
which brings the results into agreement (0.24 ± 0.05 com- 
pared to our measurement of 0.14 ± 0.1). The authors 
observed no significant evolution in the pair fraction in 
the redshift interval they studied. 

In Figure [S] we show the satellite CLF for varying host 
luminosities and redshifts. For comparison, we include 
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-21.2, Busha et al. 2011 (theory) 
-21.2, Guo et al. 2011 (theory) 
-21.25±0.5, Guo et al. 2011 
Milky Way ToUerud et al. 2008 
-21.2±0.2, Liu et al. 2011 
-21.2±0.26, Strigari et al. 2011 
-21.6 ±0.5, <2 > =0.3, All Hosts 
-21.6 ±0.5, <z> =0.3, Late-type Hosts 
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Fig. 8. — Summary of measurements of the satellite population 
of Milky Way-luminosity hosts. The legend lists the mean value 
and variation in the r band ab solute mag nitude for each study. 
The theoretical prediction from [Guo et aLI {2011b') represents the 
10 and 90 percentiles of the satellite probability distribution. Low- 
mass, low-redshift hosts in this work have a mean redshift of 0.3 
and < Mr >~ —21.6 which is brighter than the present day Milky 
Way. However, accounting for passive evolution of roughly one 
magnitude from redshift one to present, these galaxies represent 
approximate Milky Way analogs at redshift 0.3. The upper x-axis 
is shifted to the left by 0.4 mags for these hosts. 

the measu r ement of intermediate luminosity hosts from 
IGuo et all ()2011aD at redshift 0.1 along with the theo- 
retical prediction from Bll for redshift 0.3 hosts . This 
figure highlights the the strong dependence of the number 
of satellites per host on the host stellar mass and lumi- 
nosity. At the same time, there is no significant redshift 
evolution in the number of satellites. This highlights 
the important of selecting constant host mass samples 
when attempting to study trends in the satellite popula- 
tion. 

As shown by the curves, theoretical predictions using the 
techniques of Bll match our inferred satellite numbers 
well for the sample of 'All' morphology hosts dX z — 0.3, 
as discussed above. 

11. DISCUSSION 

In the Introduction we discussed a large number of phys- 
ical processes that can be tested by the spatial distribu- 
tion and number of satellite galaxies. These processes 
can broadly be divided into three sets. First how satel- 
lites are affected by interactions with their host galax- 
ies, second how host galaxies are affected by satellites, 
and third how our observations fit into the framework of 
ACDM. 

In terms of the first question, we see no evidence for 
strong physical interaction between satellite galaxies and 
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ir =-22.0±0.5, <2> =0.1 (Guo et al. 2011) 

□ -21.6 ±0.5 0.3 

□ -22.6 ±0.4 0.,3 

■ -21.8 ±0.5 0.6 

■ -22.8 ±0.4 0.6 
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Fig. 9. — Comparison of the cumulative luminosity function of 
satellites of hosts with varying luminosities. Mean host absolute 
magnitudes and redshifts are listed in the legend. The gray dashed 
and solid lines are theoretical predictions for high and low mass, 
low redshift host galaxies from this work. The satellite numbers 
depend strongly on host luminosity, while there is no significant 
trend with redshift given measurement uncertainties. 



central host galaxies on the distance scales we study 
(~ 10 — 250 kpc). This is evidenced by the lack of any 
kind of trend in the steepness of the radial profile of the 
satellite number density as a function of host redshift, 
stellar mass, or morphology or as a function of satellite 
luminosity. In a future paper, we will allow the slope 
of the radial profile to vary with distance from the host 
galaxies to allow us to compare the profile with NFW. 
We see evidence that the satellite angular distribution 
follows the host mass profile. Assuming that the satel- 
lite number density follows the dark matter halo, our 
result indicates that early-type host galaxies have light 
profiles that are aligned with their dark matter halos. 
The alignment for satellites of disk galaxies is weaker 
indicating that disk galaxies may be less aligned with 
their dark matter halos. This interpretation of mor- 
phologically dependent galaxy-halo alignment is consis- 
tent with results from weak lensing ([Mandelbaum et "ahl 
l2006f l. a nd recent theoretica l modeling of disk galaxy for- 
mation ('Deas on et al.l[20Tll ). The fact that satellites of 
early-type hosts appear in elliptical distributions which 
are more flattened than the host galaxy light profile may 
indicate that the extended dark matter halo is more elon- 
gated than the galaxy light profile at the center. Several 
simulations show that baryons can have a dissipational 
effect on the central dark matter halo, causing it to be 
more isotropic than it would have been in ab sence of 
baryons (Dubin ski 1994; Kaz a ntz i dis et aT]|2010D . 

With respect to the second question, the constant radial 
distribution and number of satellites over time indicate 



a constant accretion of satellites for all host galaxies be- 
tween reds hifts 0.8 and 0.1 for satellites brighter than 
Am < 5.5. iNewman et al.l ()201lD observed no evolution 
for bright Am < 2.4 pairs between 0.4 and 2.0. The 
combination of these results indicate that the number of 
satellites per host has been in equilibrium for at least 
half the age of the universe. From this we expect that 
host galaxies have been accr eting stellar mass at a fairly 
constant rate over this time. INewman et al.l ()201lD esti- 
mate that the added stellar mass from the bright pairs 
is roughly ^ 6% per Gyr. 

The strong dependence on host mass in the normaliza- 
tion of the satellites of early- type galaxies within a fixed 
bin of Am reflects the non-linear relationship between 
host stellar mass and halo mass. The fact that late- 
type hosts do not show a significant corresponding trend 
may indicate the effects of environment. As discussed 
above, from the morphology-density relation, early-type 
hosts are more likely to be found in groups, while late- 
type hosts exist in more isolated environments. This sug- 
gests that dark matter halo mass-luminosity relation may 
also be dependent on environment, something that has 
been largely ignored in most theoretical models. The 
relationship between host mass and satellite numbers is 
fundamental to the missing satellite problem, as clus- 
ters of galaxies have satellite numbers similar to what is 
predi cted by ACDM while isolated field galaxies do not 
(Kr avtsovll20inl: IKlvDin et al.llT999h . 

Gravitational lensing is currently the only viable method 
of measuring the ha lo mass of satellite galaxi es outside 
of the Local Group. iDalal fc Kochanekl (|2002[ ) used flux 
ratio anomalies in five quadruply-lensed quasars to es- 
timate the fraction of dark matter in satellite galax- 
ies to be between 0.6 and 7 % near the lensed images, 
demonstrating the potential of this technique to test the 
smallest mass scales of ACDM at cosrnological distances. 
More recently, iVegetti et all (|2010l 120121 ) used recon- 
structions of the mass profile of lens galaxies to detect 
low mass sub-halos near the lensing galaxies. From this 
IVegetti et al.l ()2012D estimated a mass fraction of 3.31^ g 
in satellites near lensed images. Simulations find a some- 
what l ower mass fract ion in satellites near the Einstein 
radius (IXu et al.ll200l of 0.1%. 

We can use the results from this work to make an esti- 
mate of the fraction of satellites we would expect near 
massive, intermediate redshift, early-type galaxies (typ- 
ical of lenses). Assuming a typical Einstein radius of ^ 
1" , we expect roughly 5 % of the more massive, early- 
type host galaxies in our sample to host a satellite with 
Am > 7.5, indicating an average mass fraction in satel- 
lite galaxies of a few percent in this region, on average. 
As discussed in the Introduction, the most massive sub- 
halos are stripped of dark matter in the inner regions of 
the host halo. The discrepancy between simulation re- 
sults and observation may therefore indicate the impor- 
tance of baryons in the preservation of the most massive 
satellites close to the central galaxy. 

Finally, the detection of significant anisotropy in the 
satellite angular distribution has important implications 
for using fiux-ratio anomalies to detect satellites. The 
anisotropy increases the line of sight mass between a 
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lensed image and the host galaxy, effectively making the 
host galaxy more massive in a p articular direction by up 
to a factor of 6 (Zent nei|l2006[ ). In a future paper, we 
will perform a detailed study of the predicted flux ratio 
anomalies caused by luminous satellites using the results 
from this work. 

12. SUMMARY 

We employ the host light subtraction method developed 
in Nil to study the satellites of 3425 host galaxies se- 
lected from the COSMOS field. The depth of the ACS 
images allows us to measure the satellite luminosity func- 
tion more than one thousand times fainter than host 
galaxies at low redshift, while the combination of host 
light subtraction and high resolution images allows us to 
accurately detect faint sources as close as 0.3 (1.4) arc- 
seconds (kpc) to host galaxies. Using the large volume 
of the COSMOS field, we examine trends in the satellite 
radial distribution, angular distribution and luminosity 
function. We further examine how these trends vary with 
redshift, host stellar mass and host morphology. Our 
main results can be summarized as follows: 

1. The number density of satellites of all host galax- 
ies in the sample is well-described by a power-law 
P(i?) cx with 7p = -l.ltgj. There is no sig- 
nificant deviation from this as a function of host 
morphology, stellar mass or redshift. Furthermore, 
there is no evidence for variation in 7p with satellite 
luminosity. 

2. The inference results are the same when satellite 
distances are scaled by the half-width at half-max 
of the host galaxy light profile (i?h) or by i?200 
determined from the host stellar masses. 

3. Satellites of early- type hosts, follow angular distri- 
butions which show strong alignment with the ma- 
jor axis of the host light profile, and fall in an ellipse 
that is more flattened than the host light profiles. 
In contrast, the satellites of late- type hosts do not 
show strong anisotropy. 

4. The satellite cumulative luminosity function (CLF) 
is a power-law with faint end slope approximately 
iV(L) oc i-0-5±0-i. 

5. The satellite CLF of early- type galaxies shows a 
strong mass dependence in the normalization, with 
more massive early-type hosts having significantly 
more satellites than less massive early- type hosts. 
A similar trend is not apparent for late-type host 
galaxies. This likely reflects the fact that mas- 



sive early-type hosts are more likely to be found 
in groups and exist in more massive dark matter 
halos than their late-type counterparts. This high- 
lights the mass dependence of the satellite luminos- 
ity function which is fundamental to the missing 
satellite problem. 

6. There is no significant evolution evident in the 
CLFs of satellites between redshift bins of median 
0.6 to median 0.3. Satellite numbers from the red- 
shift interval 0.1-0.4 are consistent with lower red- 
shift studies (z < 0.1) of the satellites of host galax- 
ies with equivalent masses. 

7. Predictions from subhalo abundance matching 
()Busha et all 120111) broadly agree with the mea- 
sured satellite numbers when host morphology is 
not considered. However, the theoretical predic- 
tions were not created to match trends with host 
morphology. We find that such models systemat- 
ically overpredict the abundance of satellites sur- 
rounding massive late-type hosts, and underpre- 
dict the abundance around early-type hosts. As the 
observations show smaller differences for the satel- 
lite abundance around lower mass hosts (10.5 < 
log[M^/M0] < 11.0), the theoretical model tends 
to have smaller discrepancies. Additionally, agree- 
ment is better at lower redshifts {z ~ 0.3) then 
for higher redshifts (z ^ 0.6). As discussed, this is 
likely due to errors in the model predictions, due to 
uncertainties in the luminosity function and scat- 
ter in the mass-luminosity relation at these higher 
redshifts. 
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APPENDIX 

HOST GALAXY PROPERTIES AND INFERENCE RESULTS 
SOURCE EXTRACTOR PARAMETERS AND PHOTOMETRY COMPARISON 

After removing host galaxy light from cutout images, we run SExtractor on the residual image to identify remaining 
objects. In this step we use SExtractor parameters tuned to match the photometry and completeness of the 
COSMOS photometric catalog. We could not directly the parameters listed by the COSMOS team, as their photometry 
was performed on images with O.'.'OS resolution pixels while at the time of this work, the publicly available images had 
0'.'05 pixels. 

To check the effects of using different SExtractor parameters and pixel scales, we compare our completeness and 
photometry to that of the COSMOS photometric catalog in a 4 arcmin^ cutout of the COSMOS field. In this area, we 
detect 817 objects with I-band MAG_AUTO<25. Of the 817 objects we identified, 792 have matches in the COSMOS 
photometric catalog to the same depth and within half an arcsecond. In Figure[TUl we compare the MAG_AUTO output 
from SExtractor to the values in the COSMOS photometric catalog for matching objects. The mean difference in 
the MAG-AUTO estimate is (-2.53 ± 51) x 10"^. The major outliers in the MACAUTO comparison are aU in areas 
of high object density and thus most likely due to differences in deblending. 



20 



Nierenberg et al. 



TABLE 2 

Summary of early- type host galaxy properties and results from the inference using distances scaled by ii2oo 
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1 r,+0.2 
^■"-0.2 


-o.slSi 










0.81- 


-0.81- 


9+°'2 
u.y_g 2 


3.0 


331 


11 


2±0 


1 


0.6±0.1 


-22.8 ±0.4 


1 9+0.2 
^•^-0.2 


-0.91°:^ 


> 


0.8 


< 


16 


^■"-0.2 


-0 9+'^-2 
u.y_g 2 


-| -,+0.3 
-'---'--0.2 


3.5 


322 


11 


2±0 


1 


0.6±0.1 


-22.8 ±0.4 


1 7+0.3 
^■'-0.3 


-0.81- 


> 


0.7 


< 


21 


-, 9+0.3 
-'-■^-0.2 


-0 9+''-2 
u.a_g 2 


-'-•^^-0.3 


4.0 


272 


11 


2±0 


1 


0.6±0.1 


-22.8 ±0.4 


1 S+0-3 
-■-0-0.3 


-0.81- 


> 


0.8 


< 


28 


1 4+0.3 


-0 7+°-2 
"•'-0.2 


1 7+0.3 
-'-■ ' -0.3 


4.5 


170 


11 


2±0 


1 


0.6±0.1 


-22.8 ±0.4 


2 4.+°-^ 
^■^-0.5 


-0 9+'^-2 
u.y_g 2 


> 


1.0 


< 


18 


"^■"-0.4 


-0 9+°-2 
u.a_g 2 


9 O+0.5 
^■■'-0.5 


5.0 


86 


11 


2±0 


1 


0.5±0.1 


-22.8 ±0.4 


9 9+0.7 


_n 9+0-3 
'-'•^-0.3 


> 


0.8 


< 


26 


9 7+0.8 
"^•'-0.8 


_n 7+0.3 
"•'-0.3 


^■^trr 


5.5 


31 


11 


2±0 


1 


0.5±0.1 


-22.9 ±0.3 


5l? 


-1 0+°-^ 

^■^-0.3 


> 


0.7 


< 


23 


41? 


-1 0+°-^ 


511 



^ Geometric means 
b Mean 

Satellite distances scaled by 7?200 

68 % one-sided confidence interval 

Satellite distances scaled by i?h 
f Inference performed with a Gaussian prior on 7p with mean —1.1 and standard deviation 
s Blank spaces indicate no inference on the parameter as the posterior PDF was essentially uniform 
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TABLE 3 

Summary of late-type host galaxy properties and inference results 



Am 



Nh log[M7M0]= 



Mr 



7p 



7p 



10.5 < log[M*/MQ]< 11.0, 0.1<z<0.4, 0.07 < R/R200 < 1-0, 4 < R/Ri, < 45 







„ -,-,+0.2 
'^•-'-'--0.08 






2+"-2 

^•^-0.1 






^•^-0.1 


'J-^-o.i 


-1.7+?'' 


s+°-2 


n 2+°-^ 


-1.4+0-'^ 






_1 4+0-4 
-'-•^-0.5 


o.6l?,;5 


0.8«i 


_1 4+0-4 


0.81?,;^ 


1 o+°-5 


-1 1+0-4 
-'-•-'--0.5 


1 0+°''* 
^^■'^-0.4 




-1 1+"-'* 
-'•-'--0.8 


4-^-0.5 


211 


-i.oir 


2+1 

0.91J.7 

0+3 



1.5 


274 


10.7±0.5 


0.31±0.07 


-21 


6 ±0.5 


... g 






2.0 


274 


10.7±0.5 


0.31±0.07 


-21 


6 ±0.5 








2.5 


274 


10.7±0.5 


0.31±0.07 


-21 


6 ±0.5 








3.0 


274 


10.7±0.5 


0.31±0.07 


-21 


6 ±0.5 


^^•^-0.2 


1 9+0.6 
^•^-0.9 


< 1.3 


3.5 


274 


10.7±0.5 


0.31±0.07 


-21 


6 ±0.5 


n 0+0.3 


1 9+0.6 
^•^-0.8 


< 1.3 


4.0 


274 


10.7±0.5 


0.31±0.07 


-21 


6 ±0.5 




_1 4+0-4 
^•+-0.6 


< 1.3 


4.5 


262 


10.7±0.5 


0.31±0.07 


-21 


6 ±0.5 


0.8l»i 


_1 4+0-4 
-'-•^-0.5 




5.0 


228 


10.7±0.5 


0.30±0.07 


-21 


7 ±0.5 


_1 4+0-5 


< 1.3 


5.5 


158 


10.7±0.5 


0.28±0.08 


-21 


7 ±0.5 


1 1 + ^ 

-'•4-0.7 


-1.2+?-« 


< 1.3 


6.0 


94 


lO.SitO.l 


0.25±0.08 


-21 


8 ±0.5 


2l? 


_1 1+0-5 
4-4_0.9 


< 1.3 


6.5 


39 


10.7±0.1 


0.18±0.05 


-21 


6 ±0.5 








7.0 


25 


10.7±0.1 


0.16±0.04 


-21 


6 ±0.5 








7.5 


15 


10.7±0.1 


0.13±0.03 


-21 


4 ±0.4 









11.0 < log[M*/M0] < 11.5, 0.1<z<0.4, 0.07 < ii/il200 < 1-0, 4 < R/Rh < 45 



1.5 


54 


11 


2±0 


1 


0.32 ±0.06 


-22.3 ± 0.4 






2.0 


54 


11 


2±0 


1 


0.32 ±0.06 


-22.3 ± 0.4 






2.5 


54 


11 


2±0 


1 


0.32 ±0.06 


-22.3 ±0.4 






3.0 


54 


11 


2±0 


1 


0.32 ±0.06 


-22.3 ±0.4 


5+°-5 
'^•'^-0.3 


-1.6 


3.5 


54 


11 


2±0 


1 


0.32 ±0.06 


-22.3 ±0.4 


U.D_Q 4 


-1.8 


4.0 


54 


11 


2±0 


1 


0.32 ±0.06 


-22.3 ± 0.4 


"•»-0.5 


-1.8 


4.5 


54 


11 


2±0 


1 


0.32 ±0.06 


-22.3 ±0.4 


"•8-0.5 


-2.3 


5.0 


54 


11 


2±0 


1 


0.32 ±0.06 


-22.3 ±0.4 


4^'-'-0.5 


-2.5 


5.5 


54 


11 


2±0 


1 


0.32 ±0.06 


-22.3 ±0.4 


1 9+0.8 
4^^-0.5 


-2.3 


6.0 


45 


11 


2±0 


1 


0.31±0.06 


-22.4 ±0.4 


2II 


-1.9 


6.5 


14 


11 


2±0 


1 


0.24±0.04 


-22.3 ±0.5 






7.0 


9 


11 


2±0 


1 


0.23±0.04 


-22.4 ± 0.5 






7.5 


9 


11 


2±0 


1 


0.23±0.04 


-22.4 ± 0.5 







K+0.8 

S+0.7 

5+0.6 
-0.9 

3+0.7 

-+0.7 

ui+0.7 
-0.8 

a+S'S 
-0.7 




10.5 < log[M*/MQ]< 11.0, 0.4<z<0.8, 0.07 < R/R200 < 1-0, 4 < R/Rh < 45 



1.5 


857 


10.7 ±0.1 


0.6±0.1 


-21.8 ±0.5 


2+°-2 


4-'4-0.5 


<1.0 ••• 


Q 10 + 0. 08 

'^•-'-'^-o.os 


1 9+0.5 

^•^-0.6 


2+°-4 


2.0 


857 


10.7 ±0.1 


0.6±0.1 


-21.8 ±0.5 


'^•^-o.i 


-1 2+"-''' 


< 1.1 ••• 


oq+"-os 

'^-'-'^-0.07 


-'-■^-o.g 




2.5 


814 


10.7 ±0.1 


0.6±0.1 


-21.8 ±0.5 


o.2i;!:? 


_-, + 0.5 
-■--■^-0.7 


< 1.1 ••• 






0.31- 


3.0 


726 


10.7 ±0.1 


0.6±0.1 


-21.9 ±0.5 


o.4i»:? 


„1 4+0-4 
-'--^-0.5 


<1.1 ■■■ 


2+0-4 


-I 0+0.5 


o.4lg:? 


3.5 


583 


10.7 ±0.1 


0.6±0.1 


-22.0 ±0.5 


5+°-2 

U.O_g 2 


-1 6+°-^ 
4-0-0.5 


<0.7 ■•• 


0.3«i 


_i 7+0-4 

4-'-0.6 


5+°-2 

U.O_Q 2 


4.0 


372 


10.7 ±0.1 


0.5 ±0.1 


-22.1 ±0.5 


5+°-3 
U.0_() 2 


1 7+0.6 


< 1.1 ••• 




_1 Q+0-5 
4-»-0.7 


5+°-3 

U.O_Q 2 



11.0 < log[M*/M0]< 11.5, 0.4<z<0.8, 0.07 < R/R200 < 1-0, 4 < R/Ri, < 45 



1.5 


182 


11.2±0.1 


0.6±0.1 


-22.5 ±0.5 


i+°-i 

'^--'--0.09 


-2.0l?-8 


< 1.2 ■■■ 


16+°-^ 

U.iD_Q Qg 


-1.9+r 


o.2i°:i 


2.0 


182 


11.2±0.1 


0.6±0.1 


-22.5 ±0.5 




1 a+0-5 


< 1.2 ■•• 




4-O-0.5 


0.41?,:? 


2.5 


182 


11.2±0.1 


0.6±0.1 


-22.5 ±0.5 


5+°-2 
U.O_Q 2 


9 -1+0.5 
^-4-0.8 


< 0.8 • • • 


n 7+0.3 
'^-'-0.2 


-'-"-0.4 


0.61- 


3.0 


181 


11.2±0.1 


0.6±0.1 


-22.5 ±0.5 


n Q+0-3 
^^-^-0.3 


4-0-0.5 


< 0.9 • • • 


i-ol'^;| 


-'--"^-0.4 


0.81-3 


3.5 


170 


11.2±0.1 


0.6±0.1 


-22.6 ±0.4 


^•'-0.3 


1 7+0.5 


< 1.0 ••• 




_1 4+0-4 
4-^-0.6 


Q 7+0.3 
^- '-0.2 


4.0 


139 


11.2±0.1 


0.6±0.1 


-22.6 ± 0.4 


Q 7+0.5 
"J- '-0.4 


-1 6+°-^ 


< 1.1 ••• 


o.8l°:| 


_1 1+0-6 


^•'-0.3 



Geometric means 
^ Mean 

Satellite distances scaled by i?2oo 
68 % one-sided confidence interval 
" Satellite distances scaled by Rh 

f Inference performed -with a Gaussian prior on 7p with mean —1.1 and standard deviation 0.3 

s Blank spaces indicate no inference on the parameter as the posterior PDF was essentially uniform 
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TABLE 4 
S Extractor Parameters 



Parameter Value 

Large Object Mask Point Object Mask 



DETECT.MINAREA 
DETECT.THRESH 
ANALYSIS.THRESH 
DEBLENDJvTTHRESH 
DEBLEND.MINCONT 
FILTER3FAME 
BACK.TYPE 
BACK.VALUE 
SEEING.FWHM 
INTERP.MAXXLAG 
INTERP.MAXYLAG 
WEIGHT.TYPE 



Final Photometry 



10 
L8 
1.8 
64 
lE-6 

gauss_2 . 5_5x5 . conv 
MANUAL 
0.0 
0.1 
2 
2 

MAP.RMS 



3 
3 

64 
lE-6 

gauss_l . 5_3x3 . conv 



6.5 
1.7 
1.7 
64 
0.05 

gauss_2.5_5x5.conv 



Note. — Parameters not listed are SExtractor default values. Blank spaces 
indicate that the parameter is the same for all rows. 

A full list of our SExtractor parameters is given in Table |4] 




17 18 19 20 21 22 23 24 25 

Is 14 

Fig. 10. — Comparison of our photometry with no host galaxy subtraction to COSMOS photometry in an 18 arcmin ^ field cutout. 



COMPLETENESS 

To estimate our sensitivity to low surface brightne ss companions t o host galaxies, we simulate dwarf galaxy-like objects 
in the ACS images with Sersic profiles given bv Ide Riicke et al.l (pOOl ) at a range of redshifts representative of the 
redshift range studied in this work. We then perform the host subtraction routine and estimate the fraction of objects 
we recover and the output magnitude relative to the input magnitude. The results of these simulations are shown in 
Figure 1111 For comparison we also show the results without host subtraction but with object detection performed 
using the same final SExtractor parameters as are used in the host-subtracted case. 

We choose the minimum radius at which we search for satellites where completeness begins to drop below 90% for 
early-type hosts. For late-type hosts, the minimum radius is chosen to ensure spiral arms are excluded from the 
analysis. Images after host subtraction are checked for every host to ensure that these minimum radii exclude spiral 
arms or extended residuals from disk structure (see Nil). 



FULL POSTERIOR PDFS 
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Fig. 11. — Results of completeness simulations. Red solid and blue dot-dashed vertical lines show the imposed minimum radii at which 
we study the satellite population for early and late- type galaxies respectively. On the left, the horizontal black line indicates that the input 
and output magnitudes are the same. On the right the horizontal line indicates 90% completeness. Host subtraction is essential both for 
photometric accuracy and for completeness at small radii. 



24 Nierenberg et al. 




1 2 3 4 -5 -3 -1 30 60 0.5 1.5 46 49 0.29 0.32 1 2 3 4 -5 -3 -1 30 60 0.5 1.5 46 49 0.29 0.32 

Ns 7, 1^1 (degrees) A e, (Nb/arcmin' ) <"» Ns 7, \^\ (degrees) A e, (Nb/arcmin' ) <"> 





1 2 3 4 -5 -3 -1 30 60 0.5 1.5 46 49 0.29 0.32 1 2 3 4 -5 -3 -1 30 60 0.5 1.5 46 49 0.29 0.32 

Ns 7j, (degrees) A Sj (Nb/arcmin^ ) Ns 7j, |^| (degrees) A Sj (Nb/arcmin^ ) 



Fig. 12. — Full posterior PDFs for Am =4.0 satellites of high mass, clockwise from top left: Low redshift early-type hosts, high redshift 
early-type hosts, high redshift late-type hosts, low redshift late-type hosts. Dark and light gray contours indicate regions of 68 and 95% 
confidence respectively. 



